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Abstract 

The analytic solutions of the spatially-dependent mass Schrodinger equation of diatomic 
molecules with the centrifugal term 1(1 + l)/r 2 for the generalized g-deformed Morse potential 
are obtained approximately by means of a parametric generalization of the Nikiforov-Uvarov (NU) 
method combined with the Pekeris approximation scheme. The energy eigenvalues and the corre- 
sponding normalized radial wave functions are calculated in closed form with a physically motivated 
choice of a reciprocal Morse-like mass function, m(r) = too/ (l — 5e~ a ( r ~ re " > ) 2 , < S < 1, where a 
and r e are the range of the potential and the equilibrium position of the nuclei. The constant mass 
case when 5 —* is also studied. The energy states for H2, LiH, HCl and CO diatomic molecules 
are calculated and compared favourably well with those obtained by using other approximation 
methods for arbitrary vibrational n and rotational I quantum numbers. 
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I. INTRODUCTION 



It is well known that a reasonable potential to describe the molecular vibrational spectra 
of diatomic and even polyatomic molecules successfully is the Morse potential [1-3]. There 
has been considerable current interest in the use and application of the Morse oscillator 
model to the investigation of local vibrational modes, with a particular emphasis on those 
highly excited vibrational levels which are accessible by laser spectroscopy [2]. An effective 
potential, which is the sum of the centrifugal potential term that depends on the angular 
momentum I and the Morse potential, has been used as a model for such interactions. It 
is referred to as the rotational potential. The radial non-relativistic and relativistic wave 
equations with this potential are exactly solvable for / = case [1]. However, for I ^ 
analytic exact solution cannot be determined, only numerical solutions are possible where 
several approximation techniques have been proposed and extensively used with varying 
degrees of accuracy and stability [4-17]. The most widely known approximation was devised 
by Pekeris [5] and used to obtain the semiclassical solutions [6-17]. This approximation 
[10,13] is based on the expansion of the centrifugal part in exponential terms with exponents 
that depend on an inter-nuclear distance parameter. This is why the Pekeris approximation 
is valid only for very small spatial variations from the inter-nuclear separation (i.e., for lower 
vibrational and rotational energy). Other methods that have also been used include the 
variational (V) method with the Pekeris approximation [6], supersymmetry (SUSY) using 
the Pekeris approximation [7], the hypervirial (HV) perturbation method with the effective 
potential and without the Pekeris approximation [8], the shifted 1/N expansion (SE) [9] 
and the modified shifted large 1/N approach (MSE) [9] which is very cumbersome to solve 
because there are several complex coefficients. There are needs of analytic solutions for these 
levels in order to be used for static and other applications sometimes using derivatives of 
this energy level functions. Other methods, which are semi- analytic, have also been used. 
These include, the Nikiforov-Uvarov (NU) [10] and the asymptotic iteration method (AIM) 
[12], using the Pekeris approximation, where the energy eigenvalues are obtained by simple 
transformation of the wave equation and the wave function is calculated iteratively. The 
exact quantization rule (EQR) [13,14] and the tridiagonal J-matrix representation (TJM) 
[15] which split original Hamiltonian into two parts as H = H + V where H is the part 
of the Hamiltonian that could be treated analytically while the remaining part, V, has 
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to be treated numerically. Recent techniques of approximation denoted as two-point quasi- 
rational approximation technique (TQA) [16], which is considered as an extension of the Pade 
procedure, have been used in energy hydrogenic's levels determination of Zeeman effects and 
Coulomb potential with harmonics and quadratic terms, was used for the hydrogen molecule, 
the results of Duff and Rabitz (DR) [17] and the Fourier grid Hamiltonian method (FGH) 
[18], etc. 

On the other hand, the studies of the non-relativistic and relativistic wave equations 
with position-dependent mass (PDM) for central physical potentials have attracted much 
attentions [18-29]. Such solutions are very useful in different fields of material science and 
condensed matter physics, such as semiconductors [30], quantum well and quantum dots [31], 
3 He clusters [32], quantum liquids [33], graded alloys and semiconductor heterostructures 
[34], etc. 

The NU method [35] and other methods have also been used to solve the Schrodinger 
equation [27,28], relativistic D-dimensional Klein-Gordon equation [29,36] and Dirac equa- 
tion [37] with PDM case. 

In the last few years, a considerable interest in the use and application of the mass 
dependence on the inter-nuclear distance in solving the non-relativistic and relativistic wave 
equations with various central potentials including the Morse potential have been revived. It 
provides us a reasonable first approximation for vibrational dynamics of diatomic molecules 
including the effects of anharmonicity and bond dissociation, it's incorrect behavior at large 
interatomic separations results in a relatively poor description of highly-excited vibrational 
levels, which are of much concern particularly for laser spectroscopy [2]. In other context, 
local model descriptions of molecular vibrations are required to understand the behavior 
of molecules for high overtone and combination levels [2,38]. The development of new 
experimental methods has unveiled a whole new range of molecular phenomena, including 
localization, intra- molecular energy redistribution and isomerization [39,40]. Even highly 
symmetrical molecules can develop localization at higher energy [2], a feature difficult to 
understand from the point of view of the traditional normal mode models of constant mass. 
In their pioneering work, Child et al have emphasized the need for introducing anharmonicity 
at the local level as a more natural starting point [2]. 

Since the realistic diatomic potentials are more accurately modelled by the perturbed 
or g-deformed Morse potentials, we believe that the present approach will produce much 
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more accurate information about the structure, dynamics and even the issue of energy 
resonances while dealing with the diatomic molecules. This is going to be the subject of 
investigation in the current study. In order to make the treatment very simple, we set 
out to employ the parametric generalization of the NU method combined with the Pekeris 
approximation scheme to solve the PDM radial Scrodinger equation for a generalized q- 
deformed Morse potential with arbitrary rotational angular quantum number. We propose 
and use a reciprocal Morse-like effective mass distribution function physically derived from 
the exact pseudo-spin symmetric solution of the Dirac equation [37]. Furthermore, the 
motivation for the choice is the nature of the field which is supposed to prevail between the 
interacting nuclei. This suitable choice may result in a solvable PDM Schrodinger equation 
for it's energy eigenvalues and their corresponding wave functions. Besides, it allows one to 
get a compact analytical expression and transforms the corresponding Morse Hamiltonian 
into the constant-mass problem. 

The present paper is organized as follows. In sect. 2, we briefly outline the basics of the 
NU method and supplement the parametric generalization of the method to exponential-type 
potentials in Appendix A. In sect. 3, we calculate the approximate analytic NU bound state 
energy eigenvalues and normalized wave functions of the PDM Schrodinger equation with 
the generalized g-deformed Morse potential for any /-state. The Pekeris analytic expansions 
are cited in Appendix B. In sect. 4, we calculate the numerical energy states for the various 
vibrational n, rotational / quantum numbers and q deformation parameter for four different 
diatomic molecules CO, LiH, H 2 and HCl in the constant mass limit (8 — > 0). Finally, sect. 
5 contains the relevant conclusions. 

II. NU METHOD 

The NU method has been used to solve the Schrodinger [28] , Klein-Gordon [29] and Dirac 
[37] wave equations for central and non-central potentials. Let us briefly outline the basic 
concepts of the method [35]. This method was proposed to solve the second-order linear 
differential equation of the hypergeometric-type: 

a 2 (z)u"(z) + <t(z)t(z)u'(z) + a(z)u(z) = 0, (1) 
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where the prime denotes the differentiation with respect to z, a(z) and a(z) are analytic 
polynomials, at most second- degree, and r(s) is a first-degree polynomial. Let us discuss the 
exact particular solution of Eq. (1) by choosing 

u(z) = y n {z)<i) n {z), (2) 

resulting in a hypergeometric type equation of the form: 

a(z)y:(z) + r(z)y' n (z) + Xy n (z)=0. (3) 

The first part y n (z) is the hypergeometric-type function whose polynomial solutions are 
given by the Rodrigues relation 

where A n is a normalization factor and p(z) is the weight function satisfying the condition 

[a(z)p(z)}' = r(z)p(z), (5) 

with 

r(z) = t(z) +2tt(z),t'(z) < 0. (6) 

Since p(z) > and a(z) > 0, the derivative of t(z) has to be negative [35] which is the main 
essential condition for any choice of particular bound-state solutions. The other part of the 
wave function is defined as a logarithmic derivative: 

4>{z) a(zY 1 ) 

where 

tt(z) = \ \a\z) - t(z)\ ± sj\ [a'(z)-r(z)] 2 -a(z) + ka(z), (8) 

with 

k = X-n'(z). (9) 

The determination of k is the key point in the calculation of ir(z), for which the discriminant 
of the square root in the last equation is set to zero. This results in the polynomial ir(z) 
which is dependent on the transformation function z(r). Also, the parameter A defined in 
Eq. (9) takes the form 

A = A n = -nr'(z) - \n(n- l)a"(z), n = 0,l,2,---. (10) 



At the end, the energy equation and consequently it's eigenvalues can be obtained by com- 
paring Eqs. (9) and (10). 

Let us now construct a parametric generalization of the NU method valid for any central 
and non-central exponential-type potential. Comparing the following generalized hypergeo- 
metric equation 

[z (1 - c 3 z)} 2 u"{z) + [z(l- c 3 z) (ci - c 2 z)\ u\z) + (-Az 2 + Bz-C) u{z) = 0, (11) 

with Eq. (1), we obtain 

r(z) = a - c 2 z, a(z) = z (1 - c 3 z) , a(z) = -Az 2 + Bz-C, (12) 

where the parameters c±, c 2 , c 3 , A, B and C are constants to be determined during the 
solution later. Thus, by following the method, we may obtain all the analytic polynomials 
and their relevant constants necessary for the solution of a radial wave equation. These 
analytic expressions are cited in Appendix A. 

III. NU SOLUTIONS FOR THE GENERALIZED ^-DEFORMED MORSE PO- 
TENTIAL 

Choosing the separated atoms limit as the zero of energy, the generalized g-deformed 
Morse potential has the following form [1,41,42]: 

Vm{t) = D e [q — exp(— ax)] 2 = V\ exp(— lax) — V 2 exp(— ax) + V 3 , (13) 

where a = ar e , x = (r — r e )/r e , V\ = D e , V 2 = 2qD e and V3 = q 2 D e where subscripts 1 
and 2 refer to the repulsive and attractive terms, respectively. The range of the deformation 
parameter q in the above potential was taken as q > by Ref. [43] and has been extended 
to -1 < g < or g > or even complex by Ref. [44]. The above potential contains 
three adjustable positive parameters; the parameter D e corresponds to the depth of the 
potential well, a is related with the range of the potential and r e is the equilibrium position 
of the nuclei. At r = r e , it has a minimum value at Vm^) = D e (q — l) 2 and approaches 
qD e exponentially for large r. If 1/a is somewhat smaller than the equilibrium distance 
r e , it becomes large (but not infinite) as r — > 0. Here, we refer to the above potential 
as that applicable to the diatomic molecules. For the case of multiatomic molecules, the 
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generalization of Morse potentials can be expressed by the Empirical Valence Bond (EVB) 
approach [45]. 

We now study the Schrodinger equation with the PDM for the potential form (13), which 
can be expressed as [20,26,37] 

V (-V^(r)) + 2 [E nl - V(r)\ ^(r) = 0, (14) 



m 

where m = m(r) and E nt are the real mass function and the energy eigenvalues, respectively. 
For spherical symmetry, the wave function could be separated to the following form [46-49] 

V>(r) = ^(r)lWM). (15) 

which gives the following radial wave equation: 

d 2 m! (\ d\ 1(1 + 1) 2m 



*> + ™ [r - *J - + 1? lE «< - V ' r)] \ RM = °' <16) 

where m' = dm(r)/dr and Ri(r) are the radial wave functions. It is worth noting that for 
m! = case, the above equation reduces to the well-known equation with constant mass 
used in Refs. [49-51]. Furthermore, the above equation with the use of transformation 

Ri(r) = y/mui(r), (17) 

reduces to the Schrodinger-like equation: 

d 2 u t {r) 2mE nl 

+ V eff (r)ui(r) = . 2 ui(r), «j(0) = 0, 



dr 2 j \ h 2 

l7 m" 3 fm'V m'l 1(1 + 1) 2m^. , 

V eff(r) = -7T- + T — + + -wV{r), 18 

2m 4 \ m J m r r 2 h 2 

where the interaction potential in the above equation is taken to be the generalized q- 
deformed Morse potential (13). In addition, for the bound state solutions, i.e., real energy 
eigenvalues, the wave functions Ui(r) need to be finite near r = and r — > oo. 

In order to find an approximation [4] for the centrifugal potential term l(l + l)/r 2 and for 
1/r, we apply the Pekeris approximation [5,10]. We further make the appropriate parameter 
change r = r e (l + x) and perform the expansion around x = (r = r e ) to second-order 
in x (r/r e ) at low-excitation energy (i.e., r r e ). Following the approach in Ref. [10], we 
make the convenient expansions and obtain expansion parameters a, and 6, (i = 0, 1, 2) ) 
as demonstrated in Appendix B. Thus, with the Pekeris approximations, the semi-analytic 



solution of the Schrodinger equation for the generalized g-deformed Morse potential (/ 7^ 0) 

has become amenable. For the present potential model, the exponential choice for the 

PDM allows one to transform the corresponding Morse Hamiltonian into the constant mass 

problem. Since the field between the two interacting nuclei is a short-range field of which 

the Morse potential is an example, we are tempted to set the PDM function as a reciprocal 

Morse-like ansatz having the form [52,53] 

m . 2m 5az .. 2m 5a 2 z 6m 5 2 a 2 z 2 

m = m = o , m = ~ H j-, 19 

(1-fc) 2 (1-Szf (l-5zf (l-^) 4 

where z = exp(-ax) e (0, 1) , and < S < 1. The above mass function works well in the 

present model and it results in a solvable wave equation. However, in this choice, there is 

no loss of generality when the value of the parameter 5 is taken to be small (i.e., 8 — > 0). 

Then, Eq. (18) with the help of Eq. (19), and the expansions made in Appendix B reduces 

to, 

Ul {Z) + ~z Ul{z) + zH l-5zf { W (Enl - V3) -^ + -^ + S6 ) Z 

-{^ + ^ + P6 + Q62 ) z2 } U ^ = °> (20) 

with 

s= 2_^ + 2T ?p= 26 1 _2Ta lQ= 26 2 + t%_ 
a er a cr a a z 

where and hi [i = 0, 1, 2) are given in Appendix B and ui(z) = ui(r). Defining the following 

parameters 



e n i = -Jia + 2 ~^(Vs- E^), (3, = 1 {^V 1 + 7 a 2 J + P5 + Q5 2 , 

^ = ^(^-7^)+^, (22) 

reduces Eq. (20) into the form: 

cPuAz) 1 dui(z) 1 , ^ 9 9 N , . . , 

^ + + w^f { ~' hz + ,3 ' 2Z ~ cl,) Ul(z) = a (23) 

In comparing Eq. (23) with Eq. (1), we follow Appendix A to obtain the specific values of 
the parametric constants and hence present them in Table 1 for the currently used potential 
model. Also, with the aid of Table 1, the key polynomials given in Appendix A now take 
the following particular analytic forms: 

<Z) = Znl ~ 2 (1 + ZZnl + Z, (24) 



k = [3 2 -5(2e nl + Oe n i, (25) 
t{z) = 1 + 2e nl -5(2 + 2e nl + f) z, (26) 

with t'(z) = -5(2 + 2e nl + f ) < 0, where f = ^ 1 + + | ^ - /^J. We insert the 
values of the constants given in Table 1 into the energy equation cited in Appendix A to 
obtain 



ln(n + l)5-2{n+\)^W 1 + (3 2 



(27) 



2 JK-(n + \)5 

The above equation can be written more explicitly for the rotational bound-state energy 
eigenvalues as 

2 



Em = VH 



H 2 l(l + l)a h 2 a 2 



2m r 



8m n 



n 



(n+l)5-2(n + \) + (3 2 



, n,Z = 0,l,2,... , (28) 



and after substituting the particular values of the parameters /3 1 and /3 2 from Eq. (22), 
we finally obtain the energy eigenvalues expressed in terms of the strength parameters Vi 
(i = 1, 2, 3) as 



KH(l + l)ao h 2 a 2 



E ni = V 3 + 



2mor 2 



5m 



n{n + i )5 _ 2 („ + I) ^ + !2±Ji«) + p 5 + Q5 2 + J, - *2±^) + 55 

(29) 

where 5, P and Q are constant parameters given in Eq. (21). In particular, for the constant 
mass case, i.e., in the limit 5 — > (m — > m ), we can easily reduce the above equation to, 



£n2 — 



(30) 



or more explicitly as 



h 2 



E nl = V 3 + —1(1 + 1)1 + — 2 



ar P a*rt 



2„2 



ha 
"2/7 



V2 rjt_ 

2fi 2 2/*r§ 



-1 2 



fr 2 a 2 



— 71 



(31) 
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»7 = -^7=, K = -TV 2 !"' n max < ~ "T 7i=F - 1 J > ( 33 ) 



with fi = m\mij (mi + 7712) denotes the reduced mass for the diatomic molecule. The above 
equation represents the energy eigenvalues for the generalized g-deformed Morse potential 
[42] . The vibrational bound state energy levels for s-waves (/ = 0) are 

E n = V 3 -^[l + 2n-r ] K) 2 , (32) 

with 

i 

~ / V 1 ' ah v '" max " 2\ha 
where Vi (i — 1, 2, 3) are defined after Eq. (13). Therefore, the vibrational quantum number 
n — 0, 1, 2, • • • , n max , where n max is the number of bound states for the whole bound spectrum 
near the continuous zone. Thus, n max cannot be infinite, which is reflected in the above 
condition. Therefore, the Morse eigenfunctions for real systems do not form an infinite- 
dimensional Lie algebra. 

Let us now turn to the calculations of the corresponding wave functions in the varying 
mass case. Referring to Appendix A and Table 1, we find the explicit form of the weight 
function as 

p(z) = z 2£nl (l-8z)t, (34) 

where e n \ is given in Eq. (27). The above weight function gives the first part of the wave 
functions: 

y n (z)^Pi 2 ^\l-25z), (35) 
and hence the second part of the wave functions can be found as 

<f){z) -> z £nl {l -6z)* + #. (36) 

Hence, the unnormalized wave functions are being expressed in terms of the Jacobi polyno- 
mials as 

m{z) = Kz £nl (l - 5z)^ 1+ ^P^ £ ^\l - 25 z), < 5 < 1, (37) 

where the normalization constant is being calculated using the formulas placed in Appendix 

A as 

[ r (2e nl + 1) r (g + 2) ~ (-lMn + p) (n + 1 + 2e nl + fl p 
" [ a5 e "T(n) p\(p + 2e nl )T(p + 2e nl + ^ + 2) 

1-1/2 

x 3 F 2 (p + 2e nh -n, n + 2e nl + £ + l;p + 2e ni + ^ + 2; 1 + 2e nl ; 1) , (38) 
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with e n i and £ defined in Eq. (27) and after Eq. (26), respectively, (x) is the Pochhammer 
symbols defined as 



[x 



T (x + p) 



(39) 



Jp r(x) ' 

Thus, the total radial part of the wave functions of the generalized g-deformed potential is 



ipi(r) = J\f n - [exp [— a(r — r e )]] £nl (1 — 5 exp [—a(r — r e )\) 2 + 2? 

x P^ s ^\l - 25 exp [-a(r - r e )]), < 5 < 1, (40) 

where J\f n is defined in Eq. (38). 

The constant mass case should be treated separately. To avoid repetition, we can use our 
previous calculations to find an explicit form for the weight function as [50,51] 



p(z) = z £nl exp -2^(3^ 



which gives the Laguerre polynomials: 



y n (z) -> z 



-2e„; 



exp 



J dz r 



y n+2s nl 



exp 



where y = 2J~fi[z. The second part of the wave functions can be found as 



<f)(z) — > z £nl exp - \fWtZ 



(41) 



(42) 



(43) 



Hence, the un-normalized wave functions expressed in terms of the Laguerre polynomials 
read 

Rtir) = N n (2^^ y^ exp (-|) L^-(y), (44) 

where e n i is given in Eq. (27) and y = 2^f/3^exp [—a(r — r e )] [10]. 

To demonstrate the importance of adjusting the three potential strength parameters Vi 
(i = 1,2,3) and a for real and/or complex values in any possible numerical calculation, we 
now discuss four special cases of the Morse potential given in (13) which are of much concern 
to the readers [41,54,55]. 



A. Generalized Vibrational Morse potential 



We consider the generalized vibrational Morse potential defined by [41,54,55] 
V M (x) = V ie ~ 2ax - V 2 e~ ax + V 3 , V l = D,V 2 = 2qD, V 3 = 0, 



(45) 
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and find the vibrational bound state energy spectrum as 



E n = -a 2 E 



A 



\q-n-- 



, n = 0, 1,2, • • • ,n„ 



D 1 . , 

-^rr, n max < - (2\q - 1) , 
a z E 2 



(46) 



(47) 



where the derived quantity E 



h 2 



(eV) with the following condition on the deformation 



parameter q = ^ for the final vibrational bound-state. The above result is identical to Eq. 
(30) of Ref. [55]. Also, the wave functions are calculated as 



R n {x) = A n exp j-a ^Xq — n—^jx — \e~ ax ^ L^ M n ^ (2\e- ax ) . 



(48) 



where A n is a normalizing factor. 



B. Non-PT Symmetric and Non-Hermitian Morse Case 

Following Refs. [54,55], let us assume the potential strength parameters V\ = (Ai + iBi) 2 , 
V 2 = (2Ci + l)(Ai + iBi), V 3 = and a — 1 where A 1: Bi and C 1 are real constants and 
% = y/—l. Under appropriate changes of parameters, the potential (13) turns to become a 
complex Morse-like potential: 



V(x) = -D 



+ iDe~ 



Hence, we can get the vibrational real bound state energy spectrum given by 



E n — —Eq 



1 ^ 1 

-Dki — n 

2 2 



with 



, Tl 0, 1, 2, • • • , 7l max: 



«i = jV2~ilD, n max < - {^D^D - l) , 
which is similar to Eq. (41) of Ref. [55] and the wave functions as 

R n (x) = B^yi^- 1 ^) (2K ie - x )^ dK1 ~^ e -«iexp(-x) L 2(^^-|-") ( 2Kl e- x ) , (52) 



(49) 



(50) 



(51) 



where B n is a normalizing factor. 
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C. The First Type of PT-Symmetric and Non-Hermitian Morse case 



We consider the same strength parameters as in the previous case but a — i, the potential 
(43) becomes 



V(x) = -D 



e- 2ix + iDe 



(53) 



Following the same procedure as before, we get no real spectrum for this kind of potentials: 

2 



E n = E 



1 - 1 
-D^-n-- 



with 



^ 0, 1,2, , ft max , 

1 f^e 



(54) 



k 2 = ^V^lD, n max < - (^Dy^D - l) . (55) 
in 2 \ih / 

Also, this gives the wave functions as 

Rn(x) = C n (2K 2 )'^ K ^- n ) {2K 2 e- tx )^ BK2 ^~ n) e-^xpH^QS^-i-n) ^ 



(56) 



where C n is a normalizing factor. 



D. The Second Type of PT-Symmetric and Non-Hermitian Morse case 



We consider the strength parameters V\ = u 2 , V 2 = D, V 3 = and a 
D are real constants. The potential takes the form: 

V(x) = -uj 2 e- 2iax + De~ iax , 

we get real spectrum for this kind of potentials: 



ia where u and 



(57) 



E n = E 



i Vd 



2 u, K *- n -2 



, n , 1 , 2 , • • • , n ma x j 



with 



«3 



^■y^D, n max < I \^ — ^/2ilD - 1 ) . 



h 2 \n uo 

Also, this gives the wave functions as 

R n (x) = D n (2K 3 )~(^ Ka ~^ n ) (2rize- iax )^ K3 ~^ e ~ Kze ^- iax) 

x L n {2 " 2 ; (2n 3 e- tax ) , 
where D n is a normalizing factor. 



(58) 



(59) 



(60) 
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IV. RESULTS 



Firstly, in the constant mass limit, the numerically generated vibrational bound state 
energies of the original Morse potential (q — 1, V 3 — case) for H 2 , LiH, HCl and CO 
molecules are found to be identical to those given in Table 2 of Ref. [15] (available from 
the author upon request). These numerical computations were performed using the model 
parameters [6,10,13,14,16] shown in Table 2 with the order of the eigenvalues represented by 
n (the vibrational quantum number). It is found that NU method, in the present study, can 
generate results similar to the tridiagonal J-matrix representation [15] which is relatively 
cumbersome in solving a matrix of dimension N = 100 for H 2 , LiH and HCl molecules 
to even N = 200 for CO molecule. We have also calculated the total number of bound 
states n max = 17, 24, 29 and 83 along with the whole vibrational bound-state spectrum near 
the continuous zone for the above molecules, respectively [56]. The bound state energy for 
the last state is found to be £„ max = -1.231 x 10~ 4 , -1.303 x 10~ 3 , -1.270 x 10~ 3 and 
—5.533 x !0~ 7 eV for H 2 ,HCl, LiH and CO molecules, respectively. The relative accuracy 
can be of order 10~ 5 or less (up to five significant figures) for LiH, HCl and CO molecules 
and 10~ 4 -10 -5 for H 2 molecule (up to four-five significant figures). 

In addition, the known spectroscopic values in Table 2 are used to produce the energy 
states for selected different arbitrary values of the vibrational n and rotational I angular 
momentum as shown in Table 3. We also list analogous results obtained before by other 
methods (using the Pekeris scheme) such as the NU [10], AIM [12], variational [6], SUSY [7] 
and EQR [13,14] methods together with perturbative and variational quantum mechanical 
methods like the tridiagonal J-matrix representation (TJM) [15], the shifted 1/N expansion 
(SE) [9], the modified shifted 1/N expansion (MSE) [9], the results of Duff and Rabitz 
(DR) [17], the hyper-virial perturbation (HV) [8], the two-point quasi-rational approxima- 
tion technique (TQA) [16] and the Fourier grid Hamiltonian method (FGH) [18]. The 
quality of the results is reassuring with the agreement between our numerical results and 
those generated by other methods reaching up to four-five significant digits. The current 
NU approximations to the ro- vibrational energy bound-states are slightly improved from the 
previous NU approximations probably for two simple reasons. The present calculations are 
given to four significant digits and if they rounded off to three-significant figures would co- 
incide with those given before in Ref. [10]. Throughout this numerical study, the parameter 
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conversions in Table 1 into energy units (eV) might be the second reason. The conversions 
used are 1 amu = 931.502 MeV/c 2 , 1cm" 1 = 1.23985 x 10~ 4 eU and he = 1973.29 eVx A° 
(cf. pp. 791 in [57]). A first look at Table 3 shows that the present approximations give 
results of identical accuracy like the other well-known semi-analytic variational, EQR, AIM 
and SUSY methods using the same Pekeris approximation for those fairly small rotational 
quantum numbers I. This is simply due to the Pekeris approximation [5] where the cen- 
trifugal term is being approximated to second-order in r/r e (i.e., at lower rotational energy, 
where r r e ). This is why the Pekeris approximation is valid only for very small spatial 
variations from the inter-nuclear separation. The method loses its accuracy for these higher 
ro-vibrational states. Nonetheless, our calculated vibrational energy states are in higher 
agreement with the recently calculated energy states (cf. Table 2 in Ref. [15]) for lower- 
and higher-excitation energy states. Overmore, the ro-vibrational energy states generated 
by non-perturbative NU, EQR, SUSY, AIM and variational methods are nearly same and 
slightly different from those given by the cumbersome perturbative and variational methods 
like HV, TJM, SE, MSE, TQA, FGH and DR, which are not using Pekeris approximation, 
particularly for higher-excitation ro-vibrational levels. 

The further numerical calculations of the ro-vibrational energy states for various quantum 
numbers n and / on the CO, LiH, H 2 and HCl molecules [55] are also calculated (available 
from the author upon request) using the generalized g-deformed Morse potential (V3 = 
case). The variation of the deformation parameter q in our model will produce much more 
accurate information throughout these spectra about the structure and dynamics of such 
diatomic molecules. For the case of multi-atomic molecules, the generalization of Morse 
potentials can be expressed by the Empirical Valence Bond (EVB) approach [45]. Connec- 
tions of potential functions have been extensively established for various combinations of 
pair potentials, these have been largely confined to simple potentials such as the harmonic 
[58], Leonard- Jones [59], Morse [1], Rydberg [60] and Buckingham [61] potential functions 
(cf. Ref. [62]). For example, parametric connections between the generalized Morse and 
Extended-Rydberg (ER) potential functions have been attained. Since the number of pa- 
rameters for ER exceeds those of generalized Morse by 1, therefore only Vi (i = 1,2) are 
required for converting the ER parameters into generalized Morse parameters but both the 
Morse indices are needed to obtain V\, V2 and U 3 [45,63]. The present potential functions 
have greater flexibility consisting of more parameters. There are two sets of relationship 
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between the Generalized Morse and the ER parameters: (i) for the case where V 2 < 0, and 
(ii) for the case where V 2 > 0. For example, the choice of low deformation [43] q = — 1 
provides V 2 < 0, i.e., Morse-like potential, V M {r) = D e ( e -2a(r-r e ) + 2e -a(r-r e )) with the 
position of a minimum value approaching at 14f(r e ) = 3D e and strength ratio ^ = —2. 
Whereas the high deformations q = ±5 provide the two mentioned cases for the potential 
^m(^) = D e ( e - 2a (r-re) ^ I0e -a ( r ~ re )) of minima values approaching at V M (r e ) = —9D e and 
HD e: strength ratios || = ±10, respectively. Some other numerical energy calculations can 
be estimated including the generalized g-deformed Morse potential (V 3 ^ case) [42] (avail- 
able from the author upon the request). For example, the choice of low deformations q = ±1 
provide a solution for a Morse-like potential, Vm{t) — D e (±1 — e~ a( - r ~ Te ^ 2 with a minimum 
value approaching at V^(r e ) = and 4D e , respectively. By this choice of parameters V 2 < 
and V2 > for Be-S and H-Na potentials, respectively, we can optimize the generalized 
Morse potential to ER where good correlation is seen at and near equilibrium [45,62]. This 
wide range of spectra for various deformed potential models might be necessary in fitting 
the true experimental one. This approach could be also useful to study further Ar 2 , O2, h 
and NO molecules and others. 



V. FINAL REMARKS 



To summarize, we have used a parametric generalization of the NU method derived for 
any exponential-type potential to obtain the bound state solutions of the spatially-dependent 
mass Schrodinger equation with any rotational angular momentum quantum number / for 
the generalized g-deformed Morse potential. In this paper, a suitable choice of a mass 
function has also been proposed. The present calculations include energy equation and the 
normalized wavefunctions expressed in terms of the Jacobi and Laguerre polynomials. We 
find a general equation for the ro-vibrational bound-state energy eigenvalues true for the 
currently proposed mass function given in (22) and in terms of three potential strength 
parameters Vi (where % = 1,2,3). Hence, with selected values of the parameter 5, we can 
obtain a family of solutions. The non-relativistic limit of the solution is being obtained 
by an appropriate choice of the parameter 5 — > in the mass function. Obviously, we 
can generate non-relativistic bound state solutions for the rotating Morse potential when 
the deformation parameters q — 1. The numerical application of the method to diatomic 
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molecules demonstrates that the values obtained are in high agreement with other meth- 
ods and numerical data at low-rotational excitation energy. This provides an alternative 
systematic procedure to calculate the energy eigenvalues with a reasonable accuracy and 
also considered as a suitable method in the treatment of such potentials with a varying 
mass functions as well. It is worth noting that the analytical results presented here allow 
one calculating the energy eigenvalues as well as wavefunctions in a very simple way with 
very high accuracy for lower-excitation rotational levels and for lower- and higher-excitation 
vibrational levels, probably enough for most of the applications known until now. The real 
advantages of our semi-analytic method are systematic, highly accurate, handy, easily im- 
plemented and not cumbersome as most of the other well-known methods mentioned before 
in this paper. We believe that the procedures could be easily extended to other short-range 
as well as long-range potentials. 
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APPENDIX A: PARAMETRIC GENERALIZATION OF THE NU METHOD 

We complement the theoretical formulation of the NU method in presenting the essential 
polynomials, energy equation and wave functions together with their relevant constants as 
follows. 

(i) The key polynomials: 

ir(z) = c 4 + c 5 z - [(y/cg + c 3v /c^) z - yfc%\ , (Al) 

k = - (or + 2c 3 c 8 ) - 2v^. (A2) 

t(z) = 1 - (c 2 - 2c 5 ) z-2 [(yfe + c 3v /qD z - ^ , (A3) 

t\z) = -2c 3 - 2 (yfc + c 3v /^) < 0, (A4) 

(ii) The energy equation: 
(c 2 - c 3 ) n + c 3 n 2 - (2n + 1) c 5 + (2ra + 1) (yfe + c 3v ^) + c 7 + 2c 3 c 8 + = 0. (A5) 
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(iii) The wave functions: 

p(z) =z Cw (l-c 3 z) Cl \ (A6) 

<f)(z) = z Cl2 (l-c 3 z) C13 , (A7) 

y n {z) = P^\l-2c 3 z), (A8) 

u{z) = Af n z Cl2 (l - c 3 z) Cl3 Pi Clt) ' Cll \l - 2c 3 z), (A9) 
where P^\l — 2s) are the Jacobi polynomials with 

PW\1 - 2s) = ( a + 1 >» 2Fl (_ n , 1 + a + f3 + n; a + 1; s) , (A10) 

71/. 

and N n is a normalizing factor. Also, the above wavefunctions can be expressed in terms of 
the hypergeometric function as 

u(z) = Af n z Cl2 (l - c 3 z) C13 2 Fi (-n, 1 + cio + c n + n; ci + 1; c 3 z) , (All) 

oo 

where 2 F 1 (a, 6; c; z) = f^jfk) ^ ^T^Trt^^"' ^ e can determine the normalization con- 

p=0 

oo 

stant using the condition J u 2 (r)dr = 1 and [64] 

o 

l 

J (i _ s )^-i s -i ( a> # 7 . as ) dz = 3F2 iy, a,(3 ]f i + v- 7; a) . (A12) 



(iv) The relevant constants: 

c 4 = - (1 - ci) , c 5 = - (c 2 - 2c 3 ) , c 6 = C5 + A, 

c 7 = 2c 4 c 5 - 5, c 8 = C4 + C, c 9 = c 3 (c 7 + c 3 c 8 ) + c 6 , 

2 

cio = ci + 2c 4 + 2 v /c 8 " — 1, cn = 1 - c± - 2c 4 H y'cg, 

c 3 

C12 = c 4 + v/ci, ci 3 = -c 4 + — (v^ - c 5 ) . (A13) 

c 3 
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APPENDIX B: THE PEKERIS APPROXIMATION FOR THE ROTATIONAL 
MORSE POTENTIAL 



In this appendix, we present the Pekeris approximation [5,10] performed up to the second 
order term in r/r e . We get started by making the change of parameters and coordinates as 
follows: x = (r — r e ) /r e and expanding around x = (r = r e ) to obtain: 



V rot (x) 



1 



1 + x) 



7l [1 - 2x + 3x 2 + 0(x 3 )] , 7 



1(1 + 1) 



(Bl) 



where the first few terms are sufficient for the lower-excitation rotational states since r r e . 
The corresponding rotational term expressed in the exponential form up to the second order 
is 

Vrot(x) =l{a + a ie ~ az + a 2 e~ 2az ) , (B2) 



V rot {x) = 7 



a + ai 1 — ax + 



(ax) x 



- 0(x 3 ) ) + a 2 1 - 2ax + 



(2ax) 



-0{x 3 ) 



= 7 



S y^ j a>i — a (ai + 2a 2 ) x + a 2 + 2a 2 j x 2 — 0(x 3 ) 

.i=0 



(B3) 



where (i = 1,2,3) are the expansion coefficients. Comparing (B.l) with (B.4), we obtain 
specific values for the expansion coefficients as 



a = 1 - 



a 



1 - - ) (a 

a 



2 



2 - 



, a 2 = 1 - 



a 



a 



, a = ar e . 



(B4) 



On the other hand, we repeat similar procedures to obtain an exponential expansion for the 
term 1/r around x = (r = r e ) as follows: 



r r e (l + x) r, 
or equivalently in the exponential expansion: 



' = - [1 - x + x 2 - 0(x 3 )] , 



(B5) 



- = -(bo + he- az + b 2 e- 2az ) , 
r r e 

J2bi~a(bi + 2b 2 ) x + a 2 (j + 2b 2 \ x 2 - 0{ 



X" 



Thus, comparing (B.5) with (B.6), we obtain the following expansion coefficients as 



6o = l- 



i a 



a V 2 a 



2 

a 



1 - 



a 



,h = - 



a \ 2 a 



(B6) 



(B7) 
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TABLE I: The specific values for the parametric constants necessary for the present potential. 



Constant Value Constant Value 





1 


c 2 


5 


c 3 


5 


c 4 





^0 


-¥ 


C6 


4 V 1 ^r J \j 






C8 


£ nl 


eg 


we 


ClO 




Cll 


e 


C12 


£nl 


Cl3 




01 


B 


02 


C 


£ nl 



TABLE II: Model parameters for CO, LiH, H2 and HCl diatomic molecules in our study as 

obtained from the cite d sources. 

Parameters CO [13] LiH [13] H 2 [16] HCl [13] 

D (cm^ 1 ) 90540 20287 38266 37255 

a (A )- 1 2.2994 1.1280 1.9426 1.8677 

r 1.1283 1.5956 0.7416 1.2746 

/i (amu) 6.8606719 0.8801221 0.50391 0.9801045 
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TABLE III: Bound-state energy eigenvalues (—E n i) for the H2, LiH, CO and HCl molecules (in 
eV) for different values of the vibrational n and rotational / quantum numbers in the usual Morse 
potential (q = 1, V3 = 0). 

n I This work SUSY [7] AIM [12] HV [8] DR [17] MSE [9] TJM [15] SE [9] V [6] TQA [16] 



H 2 



n n 


^±.^± 1 UUl 






A /L7RD1 

t±.^± 1 UUl 


A A7(\9 


A zL7fif) 

^±.t± 1 Uu 


A 47fifl1 ^1 

t±.t± 1 UUlOl 


A A7AQ 


A zL7^8 A 47finn8zl 




4.25880 


4.25880 


4.25880 


4.25901 


4.2592 


4.2590 


4.2590180 


4.2589 


4.2563 4.2590038 


10 


3.72194 


3.72193 


3.72193 


3.72473 


3.7251 


3.7247 


3.7247471 


3.7247 


3.7187 3.7247181 


5 


2.22052 


2.22051 


2.22052 


2.22051 


2.218 


2.2205 


2.2205369 






5 


2.04355 


2.04353 


2.04355 


2.05285 


2.050 


2.0430 


2.0528808 






10 


1.60391 


1.60389 


1.60391 


1.65265 


1.650 


1.6535 


1.6526902 






7 


1.53744 


1.53743 


1.53744 


1.53743 




1.5374 


1.5374552 






5 


1.37656 


1.37654 


1.37656 


1.39263 




1.3932 


1.3926614 






10 


0.97581 


0.97578 


0.97581 


1.05265 




1.0552 


1.0526836 






n I 


This work NU [10] 


EQR [13] SUSY [7] AIM [12] MSE [9] TJM [15] 


FGH [18] SE [9] V [6] 


LiH 





2.42886 


2.4287 


2.42886 


2.42886 


2.4289 


2.4280 


2.4288627 


2.42886 


2.4278 2.4291 


5 


2.40133 


2.4012 


2.40133 


2.40133 


2.4013 


2.4000 


2.4013352 


2.40133 


2.3999 2.4014 


10 


2.32884 


2.3287 


2.32883 


2.32883 


2.3288 


2.3261 


2.3288530 


2.32885 


2.3261 2.3287 


5 


1.64771 


1.6476 


1.64772 


1.64772 


1.6477 


1.6402 


1.6477149 


1.64772 


1.6242 
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1.62377 


1.6236 


1.62377 


1.62377 


1.6238 


1.6160 


1.6239497 


1.62395 


1.6074 


10 


1.56074 


1.5606 


1.56074 


1.56074 


1.5607 


1.5525 


1.5615114 


1.56152 


1.5479 


7 


1.37756 


1.3774 


1.37757 


1.37756 


1.3776 


1.3862 


1.3775588 


1.37756 


1.3424 


5 


1.35505 


1.3549 


1.35505 


1.35505 


1.3550 


1.3456 


1.3553770 


1.35538 


1.3309 


10 


1.29580 


1.2956 


1.29581 


1.29580 


1.2958 


1.2865 


1.2971612 


1.29715 


1.2781 
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TABLE IV: Continue. 



n I This work SUSY [7] AIM [12] NU [10] EQR [13] MSE [9] TJM [15] SE [9] FGH [18] V [6] 



CO 






11.0915 


11.0915 


11.0915 


11.091 


11.0915 


11.092 


11.0915353 11.091 


11.0915 


11.093 


5 


11.0844 


11.0844 


11.0845 


11.084 


11.0844 


11.084 


11.0843875 11.084 


11.0844 


11.084 


10 


11.0653 


11.0653 


11.0653 


11.065 


11.0653 


11.065 


11.0653334 11.065 


11.0653 


11.0653 


5 


9.79518 


9.79519 


9.7952 


9.795 


9.79519 


9.795 


9.7951838 


9.788 


9.79519 




5 


9.78833 


9.78834 


9.7883 


9.788 


9.78835 


9.788 


9.7883443 


9.782 


9.78835 




10 


9.77009 


9.77010 


9.7701 


9.769 


9.77011 


9.770 


9.7701124 


9.765 


9.77011 




7 


9.29918 


9.29920 


9.2992 


9.299 


9.29920 


9.299 


9.2991935 


9.286 


9.29920 




5 


9.29246 


9.29248 


9.2925 


9.292 


9.29248 


9.292 


9.2924786 


9.281 


9.29248 




10 


9.27455 


9.27458 


9.2745 


9.274 


9.27457 


9.274 


9.2745791 


9.265 


9.27458 




n I 


This work V [6] 


EQR [13] SUSY [7] AIM [12] MSE [9] TJM [15] 


FGH [18] SE [9] 




HCl 





4.43556 


4.4360 


4.43556 


4.43556 


4.4356 


4.4355 


4.4355522 


4.43556 


4.4352 




5 


4.39682 


4.3971 


4.39681 


4.39681 


4.3968 


4.3968 


4.3968066 


4.39682 


4.3967 




10 


4.29408 


4.2940 


4.28407 


4.28408 


4.2841 


4.2940 


4.2940628 


4.28409 


4.2939 




5 


2.80506 




2.80507 


2.80508 


2.8051 


2.8046 


2.8049687 


2.80508 






5 


2.77209 




2.77210 


2.77211 


2.7721 


2.7718 


2.7721880 


2.77230 






10 


2.68471 




2.68472 


2.68473 


2.6847 


2.6850 


2.6853673 


2.68549 






7 


2.25701 




2.25702 


2.25703 


2.2570 


2.2565 


2.2568924 


2.25703 






5 


2.22634 




2.22636 


2.22636 


2.2263 


2.2262 


2.2265969 


2.22673 






10 


2.14511 




2.14512 


2.14513 


2.1451 


2.1461 


2.1464148 


2.14656 
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